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^ ■ Context. 

I Aims. We investigate the astrometric performance of the FORS 1 and F0RS2 cameras of the VLT at long time scales with emphasis 

on systematic errors which normally prevent attainning a precision better than 1 mas. 
OO ■ Methods. The study is based on multi-epoch time series of observations of a single sky region imaged with a time spacing of 2-6 years 

at FORSl and 1-5 months at F0RS2. Images were processed with a technique that reduces atmospheric image motion, geometric 
distortions, and takes into account relative displacement of reference stars in time. 

Results. We performed a detailed analysis of a random error of positions that was shown to be dominated by the uncertainty of the 
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star photocenter determination. The component of the random error corresponding to image motion was found to be caused primarily 
^ ' by optical aberrations and variations of atmospheric PSF size but not by the effect of atmospheric image motion. Comparison of 
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observed and model annual/monthly epoch average positions yielded estimates of systematic errors for which temporal properties and 
distribution in the CCD plane are given. At frame center, the systematic component is about 25 yuas. Systematic errors are shown to 
Q ■ be caused mainly by a combined effect of the image asymmetry and seeing variations which therefore should be strongly limited to 

^ ' avoid generating random and systematic errors. For a series of 30 images, we demonstrated presicion of about 50 fI^s stable on daily, 

lyj I monthly, and annual time scales. Small systematic errors and a Gaussian distribution of positional residuals at any time scale indicate 

, that the astrometric accuracy of the VLT is comparable to the precision. Relative proper motion and trigonometric parallaxes of stars 

in the center of the test field were derived with a precision of 20 yuas yr"' and 40 //as for 17-19 mag stars. Therefore, distances at 1 
kpc could be determinable at a 4% precision if suitably distant reference objects are in the field. 

Conclusions. We prove that the VLT with FORS 1/2 cameras are not subject to significant systematic errors at time scales from a 
few hours to a few years providing that observations are obtained in narrow seeing limits. The astrometric performance of the VLT 
imaging cameras meets requirements necessary for many astrophysical applications, in particular, exoplanet studies and determination 
of relative trigonometric distances by ensuring a high accuracy of observations, at least 50 //as attained for image series of 0.5 hour. 

Key words, astrometry - atmospheric effects - instrumentation: high angular resolution - planetary systems 
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0^ 1. Introduction The best future prospects for high-precision astrometry can 

' be expected from space telescopes, but the only mission cur- 

• The availability of astrometric measurements of proper motion j-gntly planned is GAIA (Ferryman et al. |200TT) . Achieving a 

.1^ ■ and parallactic displacements at 10-100 microarsec precision single-measurement precision below 10 jt/as on V < 13 stars, it 

V provide a base for many astrophysical applications, e.g. deter- ^ill offer the opportunity to discover and study several thousands 

5^ mination of the distances to stars and their absolute luminosity, of planets (Casertano et al. 12008 ). However, GAIA operates as 

cq ^ detection of planets, microlensing studies of the mass distribu- ^n all-sky survey and cannot be pointed to a specific target, and 

tion in the Galaxy, dynamics of the Galaxy Center stars, etc . accur acy degrades rapidly towards fainter stars (Lindegren et 

The above studies imply use of very high precision astrom- aI. I2007i l. 

etry, requiring both reduction methods that are fairly insensitive in contrast, pointing to selected targets is possible with 

to major noise sources as well as telescopes fulfilling the pre- ground-based telescopes, effectively measuring distances in bi- 

cision requirements and temporal stability. The availability of nary star systems by means of optical interferometry. These 

suitable instruments however hardly matches current demand, achieve accuracies of the order of those of GAIA at V > 15. 

and is in large disproportion to envisioned future endeavors. In VLTI/PRIMA is able to measure distances between stars sepa- 

particular, this hinders programmes studying exoplanet popula- rated by 10" with 10 fias precision (Delplancke et al. 2000) with 

tions by means of astrometry, which would powerfully extend a 30 minute integration time; at a similar 20 fias h ' precision, 

and complement efforts based on other techniques and provide star separations can be measured with the Keck Interferometer 

an efficient pathway towards identifying habitable planets. (Boden et al. 1999 ). Moreover, at 30"separation in pairs of bright 

stars, an astrometric precision of 100 /ias has been achieved 

Send offprint requests to: PFLazorenko by Lane et al. (|2000|l with the Palomar Testbed Interferometer 

* Based on observations made with ESO telescopes at the La Silla (PTI), and Lane & Muterspagh (12004D and Muterspagh et al. 

Paranal Observatory under programme ID 078.C-0074 and observa- 12006D reported an accuracy of 20 ^as stable over a few nights, 

tions made with the European Southern Observatory telescopes ob- The availability of and access to ground-based high-precision 

tained from the ESO/ST-ECF Science Archive Facility astrometry facilities is however extremely limited, so that no ex- 

** Royal Society University Research Fellow oplanet detection programme has yet been established. 
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Large ground-based monopupil telescopes that operates in 
imaging mode also can significantly contribute to the detection 
and measurement planetary systems. Unlike infrared interferom- 
eters, these telescopes measure the position of a target either 
with reference to a single star or to a grid of reference stars. For 
a long time, however, astrometric measurements with ground- 
based imaging telescopes were believed to be limited by about 
1 mas precision due to atmospheric image motion (Lindegren 
I1980I I. This limitation, however, is not fundamental and rather 
reflects the performance of the conventional astrometric tech- 
nique. The first high precision observations well below 1 mas 
were obtained by Pravdo & Shaklan ( 1996) at the Hale telescope 
with a D = 5 m aperture. In the field of 90", they demonstrated 
precision of 150 //as h~'. This precision was further improved 
by Cameron et al. ('2008') with the use of adaptive optics. They 
reached a precision of 100 fias with a 2 minutes exposure and 
showed it to be stable over 2 months. 

A detailed analysis of the process of differential mea- 
surements affected by image motion allowed Lazorenko & 
Lazorenko ( I2004I I to show that the excellent results obtained by 
Pravdo & Shaklan ( 1 19961 1 represent the actual astrometric per- 
formance of large telescopes. It was shown that angular obser- 
vations with very large monopupil telescopes are not atmosphere 
limited due to effective averaging of phase distortions over the 
aperture. For observations in very narrow fields, atmospheric im- 
age motion decreases as D^^^^ (Lazorenko 2002) reaching below 
other error components. Also, the image motion spectrum can be 
further filtered in the process of the reduction by using reference 
field stars as a specific filter. Astrometric precision greatly ben- 
efits from the use of large D > 8 m apertures. 

Besides atmospheric image motion, one can list a number 
of other systematic and random error sources that could prevent 
us to reach a lOOjuas level of the precision. Many sources of er- 
ror depend on the telescope and cause long-term astrometric in- 
stability of results. To ascertain the practical feasibility of this 
new astrometric method, we have chosen the high performance 
FORS 1/2 cameras set at the VLT with excellent seeing. We have 
undertaken several tests of various time scales, ranging from a 
few hours to several years. 

The first test (Lazorenko |2006l l was based on a single four- 
hour series of FORS2 images in Galactic Bulge obtained by 
Moutou et al. ( I2003I I. It proved the validity of the basic concept 
of the new astrometric method and an astrometric precision of 
300 juas with a 17 s exposure was reached. 

In the second test (Lazorenko et al. I2007I I, we investigated 
the astrometric precision of the FORS 1 camera over time scales 
of a few days. For this study we used the two-epoch (2000 and 
2002) image series (Motch et al. 2003 ), each epoch represented 
by four consecutive nights. We reached a positional precision 
of cr = 200 - 300 yuas and detected no instrumental systematic 
errors above 30 //as at the time scales considered. The precision 
of a series of n images was shown to improve as cr/ ^/n at least to 
n - 30, which corresponds to a 40-50 //as astrometric precision. 

This paper concludes our study of the VLT camera astro- 
metric performance and extends our previous short time scale 
results to intervals of 1-5 months and 2-6 years. This covers all 
time scales required for typical microlensing, exoplanet search 
applications, and Galaxy kinematics studies. In Sect.2 we out- 
line the strategy of this study, observations, and the computa- 
tion of the star image photocenters. Astrometric reduction based 
on the reduction model (Sect.3) is described in Sect.4. The ran- 
dom errors of single measurements are analyzed in Sect.5, where 
we extract the image motion component, which proved to be of 
instrumental origin. Systematic errors in epoch monthly/annual 



average positions, and their spatial and temporal properties are 
considered in Sect.6. Astrometric precision in terms of the Allan 
deviation is discussed in Sect.7. 

2. Observation strategy and computation of 
photocenters 

As a test star field, the best choice is the field near the 
neutron star RX J0720.4-3125 whose FORSl/UTl images of 
3.3x3.3'angular size obtained in Dec 2000 and Dec 2002 by 
Motch et al. (2003) were already used in our previous study 
(Lazorenko et al. 2007). Its uniqueness is that it has the best 
history of observations available in the ESO/ST ECF Archives 
suitable for precision astrometry. Data are represented by 65 im- 
ages obtained with the B filter and obtained with a 2 year epoch 
difference, which allows for a reduction with no bias due to 
parallax. Also, the field is densely populated, containing about 
200 stars with a high light signal. We repeatedly observed it 
in Dec 2006, at integer differences of years, with FORS 1/UT2 
(lpx= 0.10"scale) and the same B filter, thus comparing model 
predicted and observed positions at three annual epochs, veri- 
fying the very long-term astrometric stability of the VLT at 2-6 
years, and computing relative proper motions used later on for 
the reduction of F0RS2 data. Observations were performed with 
the LADC optical system (Avila et al. 1997) that improves image 
quality by compensating for the differential chromatic refraction 
(DCR) of the atmosphere. 

The same test star field was imaged five times at FORS2 
(lpx= 0.126"scale) with the /^special filter with a 7" = 70 s expo- 
sure to keep star fluxes at approximately the same level as in the 
FORSl images. A one month spacing between time series was 
chosen to match the typical sampling for the observation of as- 
trometric microlensing or of the astrometric shift of stars caused 
by an orbiting planet. The availability of FORS 1 images gives 
us an opportunity to correct the measured F0RS2 positions for 
highly accurate relative proper motions determined at six year 
time intervals. This correction is critically important since elim- 
ination of proper motion from FORS2 positions decreases the 
number of model parameters, thus greatly improving the relia- 
bility of the subsequent statistical analysis. After elimination of 
proper motion obtained as shown, the positions of 5 series are 
reduced to a common standard frame with a model that fits star 
motion by parallax. Residuals of star positions (measured mi- 
nus model) are then analyzed to detect systematic errors and any 
correlations with time or magnitude. A summary of observation 
data is given in Table[T] Note the large variations of seeing which 
does not favor high precision astrometry (Sect.6). 

The primary goal of this study is the investigation of random 
and systematic positional errors of the FORS cameras. This task 
requires a careful reduction of observations, including determi- 
nation of proper motion and parallaxes with the combined use 
of images obtained with two cameras. We aim to demonstrate 
that a 300//as single measurement precision of narrow-field as- 
trometry translates to about 50//as precision for a series of 30 
measurements. 

Raw images were calibrated (debiased and flat-fielded) using 
calibration master files. Star images having even a single satu- 
rated pixel were marked and rejected for a loss, even small, of 
positional information. Positions of star photocenters X, Y were 
computed with the profile fitting technique based on the 12- 
parameter model specific for the VLT images (Lazorenko l2006l l. 
By careful examination, we developed a three component model 
that fits observed profiles to the photon noise limit. The dom- 
inant model component that approximates the core of the PSF 
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Table 1. Summary of the test field observations 



No. of nights 


camera 


spectr. 


T 


seeing, 


Date 


/images 


/unit 


band 


sec 


arcsec 


Dec 2000 


4/40 


FORSl/UTl 


B 


620 


0.49-0.78 


Dec 2002 


4/25 


FORSl/UTl 


B 


620 


0.46-0.83 


Dec 2006 


1/5 


FORS1/UT2 


B 


560 


0.55-0.65 


Nov 2006 


1/27 


F0RS2/UT1 


R 


70 


0.48-0.62 


Dec 2006 


1/37 


F0RS2/UT1 


R 


70 


0.34-0.73 


Jan 2007 


1/27 


F0RS2/UT1 


R 


70 


0.56-0.89 


Feb 2007 


1/27 


F0RS2/UT1 


R 


70 


0.45-0.63 


Mar 2007 


1/27 


F0RS2/UT1 


R 


70 


0.57-0.76 



Standard limits for seeing*: 0.47-0.78 



* Images with seeing out of these limits are affected by large ran- 
dom and systematic errors in positions (Sect.6). 
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Fig. 1. Ratio a-c/o-i„it of the mean star image size in the selected 
star sample with good images to its value in the initial sample as 
a function of flux for FORS 1 (filled circles) and FORS2 (open 
circles). 



(point spread function) is a relatively compact Gaussian with 
width parameters cr^, cr^ along x,y axes and a flux /q contain- 
ing 2/3 of the total star flux /. Two auxiliary Gaussians, each 
one multiplied by a factor or y^, are co-centered at the domi- 
nant component and approximate wings of the PSF. The model 
also takes into account the high-frequency oscillating feature of 
the PSF, computing it as a systematic discrepancy between the 
model and observed star profiles. Deviations between the model 
and observed pixel counts were found to be at the ~ 1 level. 
Determination of star photocenters is a very important element 
of the process because, as we show later on, most of the random 
and systematic errors occur at this phase. 

The precision s of the star photocenter was estimated by 
numerical simulation of random images yielding an expression 
similar to that derived by Irwin (|l9851 l for a single Gaussian pro- 
file 

FWHM I 87rcr?,/b 

This equation is valid in a much wider range of fluxes, seeing, 
and background signal 4 as compared to our former expression 
(Lazorenko et al. 20071 1. Here fluxes are given in electrons, crc 
is expressed in pixels, (p\ = 0.820 and (p2 - ^ + 0.15(o-g - 1.5)^ 
are empiric factors, and FWHM - S.IOctq is a relation valid for 
the VLT images. Due to the complex star profiles, coefficients 
(pi and 02 are not units, in which case Eq.([TJ transforms to the 
expression given by Irwin (I 1985I I. 

Computed model parameters were analyzed to detect and re- 
ject non-standardly shaped images indicating computation errors 
and actual image defects caused by blending, cosmic rays, etc. 
All images with model parameters and exceeding some de- 
liberately set thresholds were discarded. Thresholds were cho- 
sen so that the frequency of rejections was about 1 % for bright 
images. At this fixed threshold, the number of rejections grad- 
ually increased with magnitude, reaching 10-25% for faint im- 
ages, which are more sensitive to the background irregularities. 
In contrast, filtration based on caused excessive rejection of 
the brightest images, since the accuracy of the model profile is 
insufficient at high light signals and becomes comparable to the 
statistical fluctuations of counts. This gives rise to ax^ with sub- 
sequent false rejection of measurements. 

Most often, discarded faint star images have excessive size. 
This produced a selective effect seen as a systematic dependence 
of image parameters on flux /. A difference between the mean 
image size parameter era in the filted star sample and its mean 



value (Tinit in the initial star sample is typical. The systematic de- 
pendence of the ratio crG/cTmit on flux is shown in FiglT] While 
no difference is seen for bright images, at the faint end the size 
(Tg of stars selected for further processing is systematically 10- 
30% smaller. For FORSl images the effect is stronger due to 
the larger pixel scale (lower signal to noise ratio) and many cos- 
mic rays occured over the long integration time. The selection 
described here induces a similar dependence of the centroiding 
error on flux (Sect. 5). 

Further astrometric reduction revealed that some stars show 
a significant correlation between model minus observed resid- 
uals of positions and variations of seeing. This effect, detected 
primarily for relatively close star pairs, is due to the asymmetry 
of images caused by the light from the nearby star (see Sect.6). 
About 1 % of measurements subject to this effect were rejected. 

3. Astrometric reduction model 

cij(R) is a sample of / - 0,1,2 .. .N stars imaged m - 0, . . .M 
fold in the sky area of angular radius R centered on a target star 
which we denote with a subscript / = 0. In general, a){R) may 
represent only a portion of the complete sample of stars Q seen 
in the telescope FoV. Given the measured star centroids Xi,„, 
we derive, on each CCD image, the differential position of the 
target, its relative parallax, proper motion, and deviations y,„, 
from the model that may hide the astrometric signal (e.g. plan- 
etary signature). These quantities are not measured directly and 
are rather estimates of model parameters found in a certain ref- 
erence system set by the reduction model; therefore, they depend 
on it. The image m - sets the zero-point of positions, a grid of 
reference stars in this image defines the reference frame. For pro- 
cessing therefore we use relative CCD positions Xi,„ - X,„, - X,o, 

yim — ^im ^/O- 

The certain difficulty for astrometric processing is brought 
by the instability of the reference frame in time due to the atmo- 
spheric differential chromatism (e.g. Monet et al. 119921 Pravdo 
& Shaklan [T996l Lazorenko 120051 Lazorenko et al. |2007] l. vari- 
able geometric distortion, proper motions, etc. In our previous 
study (Lazorenko et al. 2007), we developed a model that cor- 
rectly handles this problem and ensures a solution in a uniform 
system with no distinction between target and reference stars. 
Here we propose a more general solution that allows an easy 
readjustment of the system of model parameters and of devi- 
ations Vim in a way that is optimal for a particular study. The 
model deals with atmospheric image motion (Sect. 3.1), geomet- 
ric distortions (Sect. 3. 2), and instability of the reference frame in 
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time (Sect.3.3). We emphasize that all the data derived from the 
differential reduction (proper motions, parallaxes, chromatic pa- 
rameters, etc.) are intrinsically relative (not absolute). This point 
is discussed in detail in Sect. 3. 5. 

3.1. Atmospheric image motion 

The variance of the atmospheric image motion in positions mea- 
sured in naiTow fields is given by Lindegren's (1 19801 1 expression 

6^ ~ (hR/Df\hRf'^T-^ (2) 

where h is the altitude of the atmospheric turbulent layer gen- 
erating the image motion and /? is a star configuration angular 
size. For a binary star, R is the star separation. Eq.(l2|i refers to 
the very narrow angle observations defined by condition 

hR < 0.5D, (3) 

otherwise we have the much worse pre- 
dicts a weak improvement of 6 with D thus limiting ground- 
based observations to a few milliarcsec precision at reasonable 
R and T. For a fictitious case of symmetric continuous distri- 
bution of reference stars around a target in a circle of a radius 
R, Lindegren obtained 6^ ~ (hR/OfihRf-'^T-^ with a stronger 
dependence on D. However, no way to practically implement po- 
tentially useful symmetric distributions was found, thus current 
estimates of 6 are based on Eq.© with milliarcsec limitation. 

We have shown (Lazorenko .2002. Lazorenko & Lazorenko 
|2004] lthat 

- any arbitrary discrete reference star distribution, of at least 
three stars, can be symmetrized; 

- symmetrization is always implemented by a standard plate 
reduction with a linear or more complex model; 

- 6^ ~ (hR/D)^(hRy^^^T^^, which suggests faster improve- 
ment of 6 with D in comparison to Lindegren's prediction 
for symmetric continuous distributions; 

- use of a large D > 8 m is required to meet condition Q for 
high stratospheric layers. 

Here we briefly summarize the concept of image motion re- 
duction based on the spectral description of this process. We 
have shown that the spectral power density G{q) of differential 
image motion in the domain of spatial frequencies q related to 
the turbulent layer is the product of two factors. The first factor 
F'(q) depends on D, the altitude h and properties of the atmo- 
spheric turbulent layer generating the image motion, and expo- 
sure time T. The second factor F" depends only on the way 
we define the "differential position" and on the geometry of ref- 
erence star distribution relative to the target. This factor is ex- 
panded in a series of even powers of q. Hence G{q) has a sim- 
ple modal structure G{q) = F'(q)Yj'^=i q^'''F'2_,{xo,yo, Xi,yi) with 
terms dependent on the distribution of stars on the sky only. 
Integration of the power density over q yields the variance of 
differential image motion 

2 (hR\\hRf'^^ 

d" ^ H2sF^^{xQ, yo, Xi, yd (4) 

^ ' .9-1 

that inherits the initial modal structure of G(q), and Hzs are 
modal coefficients. Of course, the actual value of 6^ is the sum 
over all turbulent layers with h ranging from a few to 25 km. The 
factor (hRjDf shows that large telescopes easily suppress even 
high stratospheric turbulence. Eq.Q however is valid only for 
observations at naiTow angles Q. In the case of the VLT, it holds 



only approximately at ^ < 0.5-1', which causes a problem of a 
lack of reference stars. With D < 8 m apertures, only low-layer 
turbulence is well reduced, resulting in a moderate suppression 
of the image motion. 

For optical interferometers, the atmospheric noise decreases 
as d^'^l^ (Shao & Colavita [1992| l and phase fluctuations are mit- 
igated by probing the difference of phase at ends of the instru- 
ment long baseline d. For monopupil telescopes, a fast decrease 
of image motion occurs in another way, by averaging the turbu- 
lent phase fluctuations over the aperture. The efficiency of phase 
averaging depends on the symmetry of the star configuration, 
which requires availability of the grid of reference stars. We 
emphasize that, unlike optical interferometers, monopupil tele- 
scopes are not adapted for precision measurement of the angular 
offset between a pair of stars due to the intrinsic asymmetry of 
this star configuration, for which 6^ follows dependence (|2|i. 

From Eq.(|4|i it follows that 5^ can be reduced by removing 
(filtering) the several first most significant modes s - 1,2... 
up to some optional kjl spectral mode. The residual variance 

obtained in this way depends on the first high k - 2s (k is 
even integer) active mode of the image motion spectrum and is 
therefore of a comparatively small magnitude. For the VLT, the 
gain in 6^ is a factor of 100 and more. The above possibility 
follows from the next considerations: the relative position of the 
target in image m along the x axis (here and farther on we omit 
similar expressions for y) is defined by the quantity 

Vom = aoiixo - Xi) = xo - ^ aoiXi, (5) 

tea} ieaj 

where prime indicates that index / = is omitted. Coefficients 
aoi meet the normalizing condition 2' oo, = 1 and are spec- 
ified below. Image motion in V is described by Eq.(|4]l with 
= IjZ^wi gwsifow - Y^'fiQifiwf where g„s are constants and 
fiw are values of functions f„ for a star /. Here /„ are basic func- 
tions numbered with an index w - 1,2... defined in Cartesian 
coordinates of the reference frame and formed with successive 

integer powers of the coordinates: \, x, y, x^, xy, y^ Thus 

fi\ = 1, fa = Xi, = yi, fi4 = xf ... for / star-. The expres- 
sion for that refers to some mode s involves only polynomi- 
als of X, y with a sum of powers equal to s. Indices w run from 
wi - s{s + l)/2 -I- 1 to W2 = (s -H l){s + 2)/2. The quadratic struc- 
ture of the expression for F!^^ implies that the s mode is zero 
when coefficients oq, meet conditions 

aoifiw = fow (6) 

for each w = wi . . .W2 basic function. Evidently, all modes up to 
s - k/2 - 1 vanish if oo, satisfy equations (|6]l for each s < k/2. 
For that reason, ao; are found as a solution of a linear system of 

N' = k(k + 2)/8 (7) 

equations Solution of this redundant system (because usu- 
ally > A^') is found with a supplementary condition 

2'"o,£; = min (8) 

set on the variance = 2 '^o/^? °f second item in (|5]l caused 
by centroiding eiTors for reference stars. For dense sky star areas 
the approximate expression 

FWHM k 

Arf* =- (9) 

2.36/; ^4 
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is valid where /' is the Hght flux per unit area coming from bright 
reference stars. 

Thus, the quantity V defined by (|5]l is free from the first 
modes of the image motion spectrum untill k/2 providing that 
aoi confirm to Eqs.( l6l8l ). The variance of V is 



(10) 



Because 5^ ~ due to (01, and A^. ~ R ^ according to Eq.©, 
the value of cr^ is minimum at 



= Arf, 



(11) 



which is reached at some optimal size R = /?opt (see Table |2]i of 
the reference frame (jj(R). 

3.2. Single plate reduction 

A standard plate reduction produces eff'ects equivalent to sym- 
metrization of the reference frame (Lazorenko & Lazorenko 
12004 ). Really, the basic equation of the plate model, in vector 
representation, is 



fc=x 



(12) 



where jcisa lx(A^+l) vector of x positions (including target), / 
is a matrix formed by vectors and c is a 1 x A^' vector of A^' 
model parameters c„ . We require that A^' takes only those dis- 
crete values which are defined by ^ for some optional k (k - 4 
and N' - 3 corresponds to the linear model, k - 6 and N' - 6 re- 
fer to the model that includes quadratic powers of x,y, etc). The 
least square solution of O is c = F'^f^Px where F = fPf 
is the normal A^' x A^' matrix and P is the diagonal matrix of 
weights P assigned differently for the target and reference stars. 
To comply with the image motion reduction procedure, we set 



Pii 
Pi- 







for target 

for reference stars 



(13) 



considering that &^ = + 6^ is the eff'ective variance of mea- 
surements. The best estimate of the vector x is x = ax where 
a = fF~^f^P is a projective (A^ -n 1) x (A^ H- 1) matrix that maps 
the matrix / to itself: 



«/=/■ 



(14) 



Residuals of the least square fit is the vector V = x — x, or 

V = x-ax (15) 



with the property 
V^Pf = 0. 



(16) 



The covariance matrix of V is B = {VV^} where curly 
brackets designate mathematical expectation. Using 
Eq.dTSll and considering that {xx^} = we find 

B = &^ — P~^a^ — aP~^ + aP~^a^ where &^ is a diago- 
nal matrix of elements Diagonal // elements of B are a 
variance cr? of the residuals V,-. For target {Pq - 0), all elements 
in column / = of a are zero. Hence ctq - + + which 
is equivalent to Eq.lfTOli where 

= (aP-'a'')oo = (fF-'f\n (17) 

defines the variance of reference frame component /c at the lo- 
cation of the target. For reference star / we come to a different 
expression 

aj^s^+d^- (aP-'h ^ej + 6r- (fF-'f h. (18) 



The last term is the noise from the reference frame and is a com- 
position of two components, noise from the star / itself with the 
variance crj, and nearby star grid noise with the variance A^. de- 
fined by Eq.JTTIi at the location of star /. Adding their inverse as 
weights, we find ifF~^f^)j} - + A^T^^. Hence 

cT]^a-]-Aldjl{Al + d-]). (19) 

Eqs.(fT5l), ( fT4l ) of the plate model coiTespond to Eqs.Q, ^ 
of the image motion filtration and ([8]l is the least square condi- 
tion. Therefore both methods are equivalent. However, we im- 
ply that the model ( fT2l i should include a sequence of all basic 
functions with no omission and at least k - A (linear plate 
solution) or above is chosen. The use of higher-order models re- 
sults in better filtration of image motion, though, as follows from 
Eq.®, it increases Arf. In the special case of A: = 2, differential 
positions are measured relative to the centroid of the reference 
stars; this should be avoided since it corrupts the symmetry of 
the reference configuration and greatly amplifies image motion 
to 6^ ~ ihR/Df^(hRf'^T-^ set by Q. 

3.3. Multi-plate reduction 

A single plate model ( fT2b is easily extended to the case of mul- 
tiple m - 1, . . . M images. For this purpose we specify, for any 
star /, a set of i = 1 , ... 5 model parameters ^,-.5 which are zero- 
points, relative proper motion fi^ and yu^', relative parallax tt, etc. 
Thus the position Xi„ of any star / in image m is modelled by 
Xim - TjI'^i fiwCmn + ^isV^m whcrc Vj„, are functions of time 
(of image number) coupled to and c„,,„ are model parameters 
c,v in image m. Using matrices similar to the above quantities, 
we reach 



fc = X 



(20) 



in a concatenated space formed by two types of basic functions, 
/ and V, related to the spatial coordinates or to time respectively. 
A similar expression is written fory, which requires introduction 
of the corresponding c*-'' matrix. 

Eq.(|20li is a set of (A^ H- 1) X M equations solvable by the 
least square fit with respect to N' x M parameters c„,,„ and 
(A^ -H 1) X 5 parameters ^,.5. Also, for each image m, we intro- 
duce an (A^ -I- 1) X (A^ -H 1) weight matrix P{m) with elements 
d-'^^ used for the reduction in x,y space and treated hke the ma- 
trix P of Sect.3.2. Another M x M diagonal matrix P(i) related 
to a star / takes into account the change in time of the residuals 
(x — fc)i,n. Diagonal elements of P(i} are equal to ctt^^ defined 
by Eq.lfTOli or Eq.(fT9]l depending on the star type. When star / 
measurements are unavailable at image m, the corresponding el- 
ements of matrices P(i} and P{m} are put to zero. 

Direct solution of system (l20b is, however, impossible due to 
the ambiguity between ^ and c. If some component s of param- 
eters ^,5 (for example, proper motion) systematically changes 
across the CCD, this change can be fitted by a polynomial and 
thus is not resolvable from a change in c„,„,. And vice versa, a 
change in time of some component w of Cn,,, produces an effect 
that is similar to a change in ^,5. Therefore Eq.(l20li is solved un- 
der S X N' conditions 



fPf = 



(21) 



where P is the N x N diagonal matrix of weights P,. For refer- 
ence stars, Pj are arbitrary (e.g. average of P„„ over all measure- 
ments) while Pq - for the target. The least square estimates c 
and ^ of Eq.(l20]i are found from the system of equations 

c = F-\m)fP(m)(x - ^v) 
^ = (jc - fc)P{i)y''N-Hi) 



(22) 
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solved iteratively using reference stars only with no contribu- 
tion from the target due to zero weights P()„, and Pq. Above, 
A^(i) = vP{i}v^ is a normal matrix for a star / in vector v space 
and F~^(m) is a matrix related to image m. For simplicity, 
in Eq.(l22b we omit expressions related to y. 

A simple, non-iterative solution of Eq.(l22]l exists, requiring 
only that each star measurement is available at each image, at 
constant flux and seeing conditions. In this case P, = const and 
therefore c = F~^{m}f^ Px and ^ = {x — ax)Pv^N~^{i) pro- 
viding that Pi - Pi is set. 

Given solution c and ^, we derive residuals 



V = x-fc-^v 



(23) 



which are orthogonal to the basic vectors / (at each m) and v 
(for each /): 



V'^P(m}f = 0; 



VP{i)v' = 



(24) 



Putting solution c derived from measurements of reference 
stars only into the second equation of (l22l i. we find parameters 
^o,i of the target. Hence, the variance of V at image m is 

^Im = 4 + ^0 + ^If - [v^ N-\0)v]„,n for target / = 
< = ^1- ^l^Vi^l + ^1) - [v^A^-HOv]™ ref. star /. ^ 

3.4. Reference frame quality 



In very dense reference frames, astrometric precision is limited 
only by errors of the photocenter determination if image mo- 
tion is well reduced. For sparsely populated reference frames, 
the noise Arf degrades the astrometric precision. This takes place 
even in our case of relatively low galactic latitude, -8°. Fig|2] 
shows the distribution of Arf computed for each star (processed 
as a target) and each FORSl image as a function of star dis- 
tance r from the frame center. The vertical scatter of dots for 
each star is rather large and is caused by variations of seeing and 
background. At frame center, Aj-f is minimum, about 150 yuas in 
average observing conditions. Although small, it is comparable 
to the centroiding error, s x 300 yuas, of the brightest targets. A 
good indicator of the ability of the reference frame to keep the 
output error of positions near to the precision of image centroid- 
ing is the quantity 



(26) 



which depends on the star location in the frame, its brightness, 
and in extreme cases of e^/A^^ ratio, it varies from zero to a unit 
value. In terms of y, the dependence of cr^ on the reference frame 
noise given by Eqs.lfTOli and ( fT9l ) is 



Po = (target) 
Pi > (ref. star) 



(27) 



o"/ = re;, 

where for simplicity we assumed S - Q. Although written for 
a single image only, the above equations emphasize the prob- 
lem caused by sparse reference frames with small y. Note that 
a decrease in crj with y for reference stars does not improve the 
signal-to-noise ratio (see Sect. 3. 5). 

Fig|2] shows distribution of y as a function of magnitude for 
each star processed as a target. The dependence has a specific de- 
cUne at the bright end, toy ^ 0.2-0.6. For bright targets y a: 0.5 
even at field center due to the limited number of reference stars. 
It follows that for small y, the resulting variance 0"^ significantly 
exceeds the centroiding error Sq since cr^ ~ 1/y. In this respect, 
y is a factor that specifies the quality of the reference frame. 
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Fig. 2. Characteristics of reference frame: error Arf as a function 
of distance r of the target from the field center (left panel) and 
the quality factor y (dependence on magnitude) (right panel) for 
each star and for each of 70 FORSl images. The vertical scat- 
ter of dots reflects variation of seeing and background in time. 
Reduction was performed in the standard way (Pq - 0, Po = 0) 
with = 10 (model with fourth powers of x,y) and optimal 
R = 1.5'. 

3.5. System of output data 

Due to the differential reduction, the computed parameters ^ and 
positional residuals V are relative. Weights P, define the sys- 
tem of parameters ^ of both target and reference stars. It follows 
from Eq.(l2ni which allows interpretation of ^ as the residual of 
the least square fit of some absolute parameters ^abs by basic 
functions /. Therefore what we measure are not ^abs but relative 
values 



^ = ^abs - a^abs 



(28) 



where a is a projective matrix defined similarly to a but with 
weights P put instead of P, and a^abs is a least-square polyno- 
mial fit of ^abs over reference stars in oj. Thus the reference frame 
o) and system weights P,- fully define rule ( |28T l of ^abs transfor- 
mation to ^ at A; (or A^')- Recall that weights P, for reference stars 
are arbitrary and thus so is the transformation ( |28] |. From general 
considerations we assume that P,- are equal to P„„ averaged over 
m, and Po = 0. 

Noting the similar structure of Eqs.dSTTi and ( l24l i (the first 
equation), we can apply the above considerations to the sys- 
tem of residuals V and find that it is defined by weights P,„,. 
Residuals V are related to some "absolute" residuals Vabs by the 
expression of (|28T l as 



V = Vabs - aV, 



(29) 



There is however an essential difference in treating Eqs.(l28]l and 
d29] l. The variance of Vabs depends primarily on the centroiding 
error while the scatter of ^abs is caused by the actual dispersion 
of star parameters (e.g. proper motions) within oj and can largely 
exceed random errors. Therefore, comparison of ^ values com- 
puted with different k, R or ai, will show a divergence dependent 
on the particular spatial distribution and the dispersion of ^abs- 
For example, the estimates of parallaxes computed in this test 
field for the same target, but with different k and Ropt, are scat- 
tered with a standard deviation of about ±100-200 fias indepen- 
dent of the target brightness, which exceeds the errors of paral- 
lax determination. This is quite normal and mirrors the change 
of reference system. Therefore, estimates of ^ obtained in differ- 
ent systems not can be merged into a single system, which will 
produce meaningless result. Unlike this, the change of reference 
system affects V very much less (by an order) which proves their 
merging into a single "system" (Sect.4.3). 
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System weights P„„ and P, essentially affect output residu- 
als y and model parameters ^, which is better analyzed from the 
point of view of signal detection, ze is a signal in x that generates 
some response t! in V . With regard to the target, the amplitude of 
z', according to Eq.dTSb. is t! - zsq since 0,0 = 0. With Eq.dZTll 
defining the variance of Vq, we find that the signal-to-noise ratio 
= z' /(To is T] - z yfy. Thus, while the measured signal z in Vq 
is independent of properties of the reference field, the signal-to- 
noise degrades at low y, primarily for bright targets. Now con- 
sider the reference star /. In this case z' - (l-aiijzsi according to 
Eq.dTSll. From Eq.dTSI) and the definition of 7 we find j - \- a-a 
hence z! - yze-, and r\ - z yjy- We conclude that the signal-to- 
noise ratio is equal for either type of stars, but the best 100% 
response z' in V is detected for targets. For reference stars, the 
signal decreases as ~ y, especially significant for bright stars. 

For some specific studies dealing with a full sample of stars 
(kinematics of open cluster members), uniformity of the system 
of output model parameters ^ is much desked. In this case, the 
best way is to process all stars as reference objects (f , 0, P, + 
0). The model solution ^ is then related to ^abs via Eq.d28Tl. With 
respect to proper motions, this transformation is ju,- - yfiahsj - 
Ij^eu'^'VA'abs,/ where y - I - a,, < 1 corresponds to the system of 
weights Pj. 

Untill now we have discussed reduction with reference to a 
star grid within a single isolated circular area a>{R) disregarding 
other stars in the FoV. In our previous study (Lazorenko et al. 
I2007I I, we considered reduction with multiple overlapping ref- 
erence subframes aj(i,R) each centered at each / star seen in 
the FoV. In this approach, the star / is processed, at first, as a 
target (P, = 0) measured with reference to its own local sub- 
frame a>(i,R). At the same time, this star is a reference object 
(Pj + 0) for adjacent subframes. The solution initially related to 
local frames K) is iteratively expanded to a reference grid Q 
(all FoV) and a single common system by applying a set of in- 
terlinking equations ( 1211 1. It can be shown that the final solution 
in Q does not depend on the size of the initial frames w(/, R). 
Residuals V of this solution in each image m meet conditions 
V^Pf = 0. These conditions correspond to Eq. (l24l) for refer- 
ence stars in Q. providing that P„„ = P,. Therefore a solution 
with overlapping reference subframes is equivalent to that in a 
single isolated area a>{i, 00) - or to a standard solution per- 
formed with all stars used as reference only. Consequently, no 
improvement in the signal to noise ratio is expected. This ver- 
sion of the reduction is useful for a low number of model pa- 
rameters (vector c is not used in the model), high uniformity of 
model parameters and residuals V, and a fast convergence of it- 
erations. However, assumption Pj^ = Pi used here means that 
Piin is constant in time, which is not always acceptable. 

For this study, we used standard reduction (Sect.3.3) com- 
puting the target position with weights Pq = 0, Po = to ensure 
the best response to systematic errors in V. 

4. Astrometric data reduction 

4.1. Reduction of F0RS1 images 

One of the FORS 1 sky images obtained in Dec 2000 at a seeing 
near its mean level was used as a reference. Photocenters were 
computed for 180 stars of B =18-24 mag in the central area. 
Reductions were performed with a standard model (Sect.3.3) 
that yields residuals V and model parameters ^ of a target star 
/ relative to the local grid of reference stars a){i,R). Due to 
the extremely small value of systematic errors, we tried to im- 
prove the statistics by accumulating data over all stars available. 



Therefore, reductions were repeated 180 times, processing each 
star / in turn, as a target (which till now was denoted by a sub- 
script / = 0). For that reason, the astrometric precision varied 
depending on the distance r of the star / from the frame center. 
For the brightest stars, this occurred due to y decreasing from 
y » 0.6 at the frame center to 0.2 at the periphery (Fig|2]i with 
a corresponding error increase by 25-50%. Computations were 
carried out with all k from 6 to 16 and all R from some R^i„ 
(Table |2]i that provides the minimum number of reference stars 
needed at a given k, to the maximum P^ax = 2.3'. The first run 
was performed with zero image motion and afterwards compu- 
tations were repeated with the actual estimate of S (Sect. 5). 

Table 2. Minimum and optimal radii of reference fields 



k 




^opt 


6 


40" 


50" 


8 


50" 


70" 


10 


60" 


90" 


12 


70" 


110" 


14 


80" 


140" 


16 


90" 


-140" 



We assumed a 10-parameter model for ^ with zero-points, 
proper motions yu', atmospheric differential chromatic pa- 
rameter p, the LADC compensating displacement p', and with 
no parallax for integer year differences between epochs. Four 
extra g^^, g„, gy^, and gyy model terms for each star were ap- 
plied to compensate for a strong, over ±200 px jittering of im- 
ages. Jittering induced a large signal in V, clearly correlated with 
telescope displacement A^, A,, along the x, y axes. This effect in- 
creased with R, reaching several milliarcseconds at P a; Pmax, 
< 8. The jittering makes the reduction difficult, and is the rea- 
son that we discarded a linear (k - 4) model reduction. 

Jittering moves the star field with reference to the cam- 
era optics. This causes a change in the optical distortion (say 
along x) at some point x,y from its initial value (^t(x,y) to 
O.V (jc, y) + (d^_^(x, y)/dx)Ax + (d<i>Jx, y) /dy)Ay. A similar expres- 
sion is valid for the distortion (!>y{x,y) along the y axis. Naming 
the partial derivatives used here as g„, we come to the expres- 
sion 

O'^ix, y) = (&^(x, y) + g„A^. + g,,Ay 

<D;(x, y) = <by(x, y) + gyAx + ^yjAy, ^ ' 

equally applicable for the correction of positions. Assuming that 
optical distortions are stable over the observing period, gt, terms 
are included in, and found as components of ^. 

An example of the measured and model star motion over the 
CCD surface is shown in Figl3]for a red (B = 20.8, B - P = 2.8) 
star with proper motion /i-* - -1.32 + 0.06 mas yr"', yu-' = 
10.74 + 0.06 mas yr"', and trigonometric parallax (all relative) 
n = 0.286 + 0.053 mas which was computed later based on 
FORS2 images (Sect.4.2). This graph is actually rather simpli- 
fied because it refers to positions corrected for polynomial fc 
and for jitter related terms. The intricate shape of the track is 
due mainly to the DCR shift of images within a single series. In 
the blue band, this motion exceeds ±10 mas while in the red fil- 
ter the effect is an order lower. This makes it clear, for instance, 
why reference star displacements should be taken into account 
when processing B images. 

The dependence on magnitude of the internal precision of 
proper motions is shown in Fig|4] These estimates correspond to 
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Fig. 3. Example of B = 20.8 star motion over a CCD traced 
for 6 years with FORSl and 5 months with FORS2: measured 
positions (open circles with eiTor bars) and model track (solid 
curves). Reduction was performed with parameters k - 10 and 
R = 1.5'. 




re 

"a 
o 



400 
350 
300 
250 
200 
150 
100 
50 




proper 


motions [uas yr"^] 




^ A 




parallaxes [|j.as] \ 1 £ 






A 






\ A 


zi? 




°» 


















O 




FdRS2 : V 




^ o _ 






FO 








RS1 - 



16 



17 



18 



19 20 
magnitude 



21 



22 



23 



Fig. 4. Astrometric precision of relative proper motions (open 
circles) determined from FORSl images with a six year time 
span and trigonometric parallaxes (triangles) derived from 
FORS2 images with corrections based on FORSl proper mo- 
tions. Reductions performed with A: = 10 and R = 1.5'. 



the formal least squares precision and take into account compo- 
nents s, 6, and Ajf of the total random error cr. Due to the large 
time span between epochs, proper motions were derived with 
high precision, reaching 20 juas yr"' for bright stars. Systematic 
errors (Sect. 6) degrade precision little since their contribution is 
small in comparison to random errors (see Sect.4.2). 



4.2. Reduction of F0RS2 images based on F0RS1 proper 
motions 

For processing of FORS2 images we used high-precision proper 
motions yu, derived from the reduction of FORS 1 images on a 
six year time base. At five-month spacing, they ensure very ac- 
curate corrections. The possibility to use proper motions that fit 
image measurements of the other camera, different wavelength, 
and at distant epochs, however, is not evident and should be 
used with care. On the other hand, the precision of FORS2 posi- 
tions obtained in this way is an indicator of the actual accuracy 
of FORSl proper motions. Most important, the elimination of 
proper motions from model parameters ^ essentially increases 
the DoF (degrees of freedom) of least square residuals in tempo- 
ral subspace thus making further study of systematic errors more 
reliable. 

The reduction was started by finding the FORS2 image ob- 
tained at normal seeing and best matching the star content of 



FORSl images. This image at epoch To was used as a reference 
for the reduction of all FORS2 images. In most cases, the dif- 
ference in star content from the two cameras occured for a small 
gap between two CCD chips of FORS2 and saturation of bright 
stars in the R filter. In this way, 152 common stars were selected 
for further reduction, which started from applying coiTections 
for proper motion occuring in star positions between epochs To 
and r,„. This was performed taking into account the singular- 
ity of astrometric reduction according to which the i-th target 
position (and proper motion yu,) is related to a particular subset 
a)(i, R) of reference stars j G a>(i, R), whose unique model pa- 
rameters (denoted by jj'-j in contrast to yu,) are valid within this 
subset only. Therefore, to conserve the reference system, a re- 
duction of FORS2 images for a target / was performed with ref- 
erence to the same frame tij(/, R) as used for FORS 1 . Also, we 
applied weights P, that are same as those involved in the reduc- 
tion of FORS 1 images, that is, using average light fluxes in the 
blue filter. Thus, each reference (with respect to target /) star 
j € aj{i, R) positions were corrected by 



-MijiTm - To). 



(31) 



This complicated procedure is due to the necessity to conserve 
the system of model parameters when processing different sets 
of images. Violation of this principle immediately destroys the 
accuracy. Thus, direct application of corrections Aym = -fj.j{Tm- 
Tq) to all measurements of reference stars in <y(/, R) is incorrect 
because these fij are related to their own frames a)(j, R) which 
differ from a){i,R). Mismatch of these areas and even a small 
inconsistency of the reference star ensemble sometimes result 
in large 1000-5000 //as epoch residuals. Even use of proper 
motions fi'^j did not ensure complete identity of Lo{i, R) related 
to each camera due to inavailability of some FORSl stars in 
FORS2 images. 

The reduction model included zero-points, chromatic param- 
eters p, p', and trigonometric parallaxes n. Formal random pre- 
cision of F0RS2 parallaxes for stars of different brightness is 
given by Fig|4] For the best stars, relative parallaxes are deter- 
mined with a precision near to 40 //as, which means that dis- 
tances at 1 kpc are measurable with a 4% precision. A few large 
upward deviations in Fig|4]for some stars are caused by a low 
number of measurements (oversaturation of bright images, or 
position in the gap between two chips of the camera) or by the 
peripheral position of stars and thus low y (large A,f ). 

The precision of parallaxes is increased by use of FORS 1 
proper motions, allowing us to exclude a component fi'' from 
model parameters ^, removing in this way a strong coiTelation 
between n and fi \ In other cases, the expected precision of par- 
allaxes from the 5-month series of observations is 200 //as only. 

Systematic errors, of course, affect the accuracy of both 
proper motion and parallax determination. In Sect. 6 we show 
that the systematic eiTor for targets near frame center is about 
25 //as, or a half of the random error of epoch average positions 
for bright stars, at months to year time scales. Translating this 
estimate to parallaxes, we find that systematic errors contribute 
approximately ±20 //as to each star parallax and ±10 /(as yr"' to 
proper motions irrespective of the star magnitude. 

4.3. Merged residuais {V) 

While processing, we computed residuals V for each star /, each 
image m, all reduction modes k from 6 to 16, and several ref- 
erence field sizes R, including /?opt- The best precision residuals 
computed atR- Ropt we denote as Vi,„{k). In practice, however. 
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at any A: and the use of (V) instead of V is of low efficiency. 
For bright stars (FiglSj?), the precision of {V} is better due to the 
averaging of the reference frame noise. 

5. Random errors 

5. 1 . Calibration of tlie image centroiding error s dependence 
on flux 

In this Section, our study was carried out with images obtained 
in a narrow seeing range of 0.47-0.78", which includes almost 
all FORSl and about 80% of FORS2 images. The use of images 
out of this range leads to a noticeable increase of random errors. 



Fig. 5. Image-to-image change of FORS 1 residuals V for a) a 
faint 22 mag star (e = 20mpx - 2000;uas) and b) bright 19 mag 
star [e = 3mpx = 300/ias); various line types refer to reduction 
parameters k from 6 to 16. Computations were made at Ropt. 



we do not require multiply defined residuals but rather a single 
set of residuals which for a particular star / is the best estimate of 
the planetary signal at the moment of image m exposure. For that 
purpose we merge Vi„,(k) into a single system, the possibility of 
which follows from the discussion in Sect. 3. 5. 

Let us consider Fig|5]that presents an example of the image- 
to-image change of Vim(k) computed with different k for two 
stars of different brightness. Residuals corresponding to different 
k are seen to be highly correlated, especially for a faint star, and 
fluctuate near their average, (y)„„ being a function of m. Recall 
that according to (|29ll, Vi,„{k) = Vabs,™ - Z'j a;/^) Kbsj-m where 
aij(k) refer to k used. Therefore (V)™ = Vabs,™ - 'Z'/a/jW^b^jn, 
where (fl,j) is an average of aij(k) with respect to k. Hence 
Vimik)-{V}im = i;}[fl,/(fe)-<fl,7)] Vabsjm- The Variance of this dif- 
ference, neglecting the second term, is Z'jlaijik)]^^^^^^, or A^(A;) 
at k given. Thus, the standard deviation of Vi„,{k) - {V)i„, de- 
pends on Arf almost linearly. This approximation is confirmed 
by actual data, as shown in Fig|6] 
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Fig. 6. Standard deviation (scatter with respect to k) of V{k)-{V) 
for each FORSl star as a function of Ajf (open circles) and a 
linear approximation (solid line). 

Given 6 sets of Vi„{k) corresponding to A: = 6... 16 for each 
target /, we merged them into the weighted average {V)i,„ using 
weights AjTf^. Along with V, merged residuals (V) were tested 
for the presence of systematic errors (Sect. 6). As explained in 
Sect. 3. 5, the merging is not applicable to model parameters ^. 

For faint stars, the relative amplitude of V fluctuations near 
(V) is insignificant (Fig|5^) since A^ <^ s. Therefore (V) s; V 
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Fig. 7. Initial (dashed lines) and corrected (solid Unes) resid- 
uals of the error expansion ( l32l i computed at minimum refer- 
ence field size /^min- Different lines correspond to reduction with 
k = 6 16. The case of A: = 6 is marked by open circles (devi- 
ating pair of lines for FORS 1 ). 

The use of stars of different brightness to investigate sys- 
tematic errors requires careful calibration of the dependence on 
flux of the image centroiding error e ([TJ. For calibration pur- 
poses, the best residuals are V computed at the minimum pos- 
sible R = Rmin (Table |2]i since they contain negligible input of 
atmospheric image motion 5 = 0. Using the variance crj of resid- 
uals V computed for each A at /? = ^min, we can find the residual 
discrepancy of the decomposition (|25] l into error components 

AL = ^ - - + [v''N-'(i)v]mm (32) 

for each target i. The change of this quantity with star magnitude 
is shown in Fig|2]by dashed lines for each A. All curves corre- 
sponding to different A modes closely follow a common depen- 
dence with little scatter. The anomalously large deviation seen 
for the FORS 1 camera at A = 6 originates from the large jit- 
tering of images which was not completely compensated by the 
reduction. For high A modes this effect is well removed. While 
for bright stars, discrepancies A^^.^ are fairly small, at the faint 
end we note a systematic negative bias caused by incorrect mod- 
elling of errors. This bias almost does not depend on A and is ap- 
proximately proportional to e^. Therefore we assumed that this 
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discrepancy is caused by insufficient validity of model ([T]) for 
e^, which requires an additive correction ifs^ with a coefficient 
ip independent of flux. A similar correction should be applied 
to also. Correction factors + tp to e and Arf computed in 
0.5 mag flux bins are shown in Fig[8]as a function of magnitude. 
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Fig. 8. Correction ^\ + ip for the model error of photocen- 
ter measurements s. Estimates for each k reduction parameter 
(dashed lines) and their average (thick solid lines). Opened and 
filled circles reproduce the ratio of image size (FiglT) in the se- 
lected and initial star sample 



A change of ^\ + with brightness in Figl8]is similar for 
both cameras. A negative trend over a 4-5 mag range of bright- 
ness reproduces the dependence of star image size ctg on flux 
(FigHJ and therefore probably is a consequence of selective fil- 
tration based on star profile parameters when star images with 
excessive size were discarded (Sect.2). The use of more compact 
images in comparison to the initial star sample, of course, results 
in an improvement of the effective centroiding error e observed 
in Fig|8] A similar improvement of precision for the brightest 
images occurs for the selection based on;^^^ criterion (Sect.2). 

Averaging with respect to k produced final estimates ^\ + ip 
shown in Fig|8]by solid lines. With these corrections, residuals 
(|32] | have been recomputed yielding new discrepancies A^^^ with 
much smaller magnitudes (Fig]?] solid lines). Having found the 
caUbration factor + tp, we can correctly estimate o-y at any 
R > Rmin simply by adding the image motion variance 6^: 



o-l = (1+ <p)s^ + (1 + v)Al + 6j + A^ - [y'N-'(i)y]„ 



(33) 



The term A^^.^ is used to take into account the dependence of 
(p on flux which originally was considered constant. This also 
compensates, at least statistically, the use of a single correction 
factor for both and A^. 

The validity of above calibration is illustrated in Fig|9]where 
we compare the astrometric precision of a single photocenter 
measurement restored from observations with its model predic- 
tion s ^\ + ip in the case of reductions with R - 1.5 'and k - \Q. 
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Fig. 9. Astrometric precision of a single photocenter measure- 
ment: observed (open circles) and model estimate + p 
(filled circles) for FORSl (large symbols) and FORS2 (small 
symbols) as a function of magnitude. 

The measured astrometric precision, for each star, was computed 
based on the observed variance cr^ of V (mean in x and y), 6^ 
derived in Sect.5.2, and representation ( l33T l. These results, as 
for model values s + p for each target, were averaged over 
all data available. Fig|9] shows a good match of the observed 
and model precision over wide range of magnitudes. This graph 
matches well our previous results for FORS 1 based on a reduc- 
tion technique with overlapping reference frames (Lazorenko et 
al. l2007l l. 

We emphasize that both e ^]\ + ip and e are estimates of the 
actual precision of the photocenter determination. The difference 
is that the first one refers to the star sample affected by selection 
while e is related to the imaginary sample of FORS images with 
no defects. In spite of the small value of p, the subsequent study 
of image motion and systematic errors greatly favours its use 
since it allows us to incorporate large amount of data from faint 
stars. 



5.2. Image motion 

Taking advantage of the availability of a well calibrated image 
centroiding error, we used Eq.([33]l to extract the image motion 
component 6 at various R. This equation was solved numerically 
for each star taking into consideration the fact that the refer- 
ence frame noise Arf is a function of e and 5. The results av- 
eraged over all stars available and computed for each k and R are 
shown in Fig [10] Comparing estimates obtained for both cam- 
eras, one may note the similarity of results in spite of the dif- 
ference in pixel scale, number and spacing of epochs, different 
reduction model parameters, different method of reduction and, 
especially, an 8-fold difference in exposure T (600 and 70 sec 
for FORSl and F0RS2 respectively). The last aspect raises a 
doubt about the validity of relating the measured image motion 
to atmospheric turbulence. 

At each fixed k, 6 estimates were fitted by a power law 



S = BR" 



(34) 



assuming that R is given in minutes of arc. Fitting parameters B 
and b are given in Table[3]for the first few k modes only since the 
results for ^ > 12 are too uncertain. Excessive estimates of B (in 
comparison to FORS2) found at ^ = 6 and k - E could be due to 
the residual effect of large image jittering of FORS 1 images. For 
comparison, the table contains Z?a and b^^ coefficients of Eq.(|34li 
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observed image motion at T > 70 s is not due to atmospheric 
turbulence since it does not decrease as r '^^ and therefore is of 
instrumental origin. Very likely, it does not depend on exposure, 
at least for T > 70 s. Due to domination over the intrinsic atmo- 
spheric image motion, the last component not can be extracted 
from the present data. In Sect. 6 we suggest that the image motion 
detected probably is caused by star image asymmetry in combi- 
nation with variations of the PSF. 



6. Systematic errors at monthly/annual epochs 
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Fig. 10. Image motion <5 as a function of reference field size 
R computed with k = 6 ... 16 (symbols of different type) and 
corresponding fits (lines) by a power law ( |34] |. 



expected for differential image motion caused by atmospheric 
turbulence. These values were obtained by scaling model coef- Ve - ^ VimcrJi^ / ^ 
ficients (Lazorenko & Lazorenko I2004I I that refer to typical at- tnee mee 

mospheric conditions at Chilean observatories, to current expo- 
sures and telescope aperture. The atmospheric-related amplitude 
Ba is much smaller than the observed one, especially for FORS 1 
with a 3-8 fold discrepancy. Such a large difference suggests that 
we are measuring an effect not related to the atmospheric turbu- 
lence. 

In a pilot study of FORS2 astrometric performance, 
Lazorenko (2006) estimated Eq.(l34li parameters B' and b' us- 
ing a single night observation series with T - 11 s exposure. 
Coefficients B' reproduced in Table [3] are approximately twice 
as large as in this study, possibly due to the different technique 
of reductions, which now takes into account DCR displacement 
of reference stars. 



In our former study (Lazorenko et al. |2007| l. systematic errors in 
positional observations with the FORS 1 camera were shown to 
be about 30 juas. The detection of such weak signals presents a 
certain difficulty and limits our search to characterization of the 
error component invariable within each monthly/annual epoch 
of observations. Systematic signs in observations may appear for 
effects not described by the reduction model. The most troubl- 
some are long-term instabilities which differently affect images 
at distant epochs, e.g. changes in VLT optical aberrations, star 
colours, actual PSF shape, variable background gradient due to 
light from nearby stars, etc. 

6.1. Epoch average residuals 

The quantities best suited to this study are the epoch average 
weighted residuals 



(35) 



computed for each star at each monthly/annual observation 
epoch e. A set of these epoch points Ve was investigated to detect 
and characterize systematic errors. The epoch normal points Ve 
are formally characterized by variances 



De 



(36) 



equal to the cumulative weight of individual residuals of images 
m € e available at epoch e. However, variances De describing 
the actual scatter of the normal points Ve are below De as a con- 
sequence of the least squares fit. Thus, assuming a normal law 
for the distribution of observation errors, from Eq.(l35]l we find 



Table 3. Coefficients of Eq.([34li: derived in this study B[yuas], b; 
predicted B^, b^ by atmospheric model (Lazorenko & Lazorenko 
I2004I I: and B', b' obtained from a single series of F0RS2 images 
(Lazorenko ISOOgj i 
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B mm' itn 



^ im ' 



(37) 
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0.86 


27 


1.8 


159 
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1.40 
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where B„m' are diagonal elements of a covariance matrix 

B{i) = p-\i) - v^N-'^{i)v (38) 

of residuals y,„, for /-th star and P{i) is a diagonal matrix with 
elements cr^^ introduced by Eq.JTOli. It follows that De - De 
only when the second item in Eq.(l38]l is zero. In practice, instead 
of Eq.d37Ti, it is convenient to use the expression 



De = OiDe 



(39) 



'■ Could be biased due to the residual effect of large image jittering. 



In all cases, the measured powers b of Eq.(l34ll are signif- 
icantly below their predicted values ba_. We conclude that the 



where 6e < 1 is a quantity numerically computed for a particular 
distribution of /-th star observations over time. Although matri- 
ces N~^{i) are unique for each star, a minor difference in Be for 
different stars often can be neglected. Of course, we have two 
sets of 6e values related to x and y axes. 

It is difficult to suppose that systematic components follow 
exactly the parallax and proper motion displacement of stars. 
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Fig. 11. Monthly normal points (average residuals) Ve for stars observed with FORS2. Symbol size refers to the expected precision 

1 /2 

of normal points D^' , which changes from better than 100 /vas (large circles, usually brightest stars), to 170 juas (middle), and 

1 /2 

260 //as (small circles). Corresponding 3cr scatter areas are shown by dashed ellipses with semiaxes 36eD^' whose size depend 
primarily on seeing (best in December 2006 and worst in January and March of 2007). No large dispersions are seen. 



Therefore, after a fit in time, systematic errors add an extra scat- 
ter to epoch residuals Ve, which is detected as an excess in the 
expected value of the variance D^. This excess we find below 
based on the well-calibrated (Sect. 5) model of the stochastic- 
dependent component of the variance. 

A good idea of analyzed epoch average residuals Ve is given 
by graphs of Fig[TT] This plot shows the typical distribution of 
monthly normal points Ve^\ V^'^ in the x, y plane for each month, 
for stars observed at FORS2, and the reduction with A: = 10 
and R - Ropt- This distribution of Ve is typical also for pro- 
cessing with other k at Ropt due to the high degree of correla- 
tion between these sets of residuals (Sect.4.3). Normal points 
are shown by open circles of three size grades which refer to the 

1 /2 

precision De of normal points better than 100 juas (largest cir- 
cles, usually brightest central stars), 170 //as (middle size), and 
260 /ias (small circles, faint or peripheral stars). Ellipses (dashed 

1 /2 

curves) with semiaxes WeDJ mark 3cr scatter limits expected 
for least square residuals. The different scatter of Ve for differ- 
ent epochs is caused primarily by seeing (Table[T]), which is best 
for the second (most compact location of dots) and worst for the 
third and last epoch. For the same reason, only a few points with 
dI^^ < 100 //as precision (large signs) are seen for the middle 
and the last epoch since only the best stars are measured well 
at bad seeing. Most normal points of each precision grade are 
inside of the corresponding 3cr limits with no wide dispersions. 

Next, we considered the frequency distribution of Ve- These 
data however are not uniform in precision for different light 
fluxes from stars. In order to exclude a dependence of the preci- 
sion on brightness, we introduced the dimensionless normalized 
quantities 

Ve = VeD-'l\ (40) 

The standard deviation of Ve, according to ( l39b . is equal to 9e 
and therefore does not depend on brightness. The frequency 
distribution of Ve was formed cumulating data of all epochs, 
model versions with k - 6. ..16 at 7? = 7?opt, for both axes, 
and using bright B < 21 stars in the central frame area r < 
1'. Histograms obtained and their Gaussian approximations are 
given in FigfT2h.b for FORSl and FORS2 respectively. These 
histograms are compared with the theoretical distribution of Ve 
in the case of zero systematic error and taking advantage of the 
fact that Ve follows a Gaussian distribution with the variance pa- 
rameter {6)- equal to 6^ averaged over the epochs. 

In the case of FORS 1 , from Eqs. (|36l - l38] l we find that typical 
values of 0^ are 0.17, 0.48, and 0.27 for epochs e = 1,2, 3 respec- 
tively, with small variations depending on observing conditions 



of the particular target. Therefore {6) - VZi 6'e/3 - 0.56 is a 
cr-width parameter for a theoretical Gaussian distribution. The 
observed distribution of Ve for FORS 1 is slightly wider, with a 
0.61 cr-parameter and a few large residuals (FigfT2h). 

In the case of FORS2, 6^ varies for different epochs between 
0.22 and 0.76, with an average {0) - 0.78. The observed distri- 
bution is much wider, with a 1.17 cr-width parameter and a sig- 
nificant widening of wings (FigfT2b). The observed distribution 
of monthly normalized residuals clearly indicates the presence 
of large systematic eiTors which is discussed later on. 

The histograms in Fig [12] are sensitive to systematic er- 
rors providing their magnitude is comparable to the preci- 
sion of epoch normal positions, about 50-200 //as. These his- 
tograms however are to be considered primarily as illustra- 
tive. Numerical characterization of systematic eiTor is found un- 
der the assumption that its value Ag for a given star / is con- 
stant within each monthly/annual observation series e. In this 
case all measured residuals y„„ in images m e e are system- 
atically biased by a constant Ae- Therefore yf(measured) = 
Ve{a.t zero systematic errors)-!- A;, and the expectation of the vari- 
ance of measured normalized residuals is Ve - Ol + A^d^D^K 
The second item describing the input of systematic errors in Ve 
dominates for bright stars. This component was computed for 
each star and averaged to derive statistically reliable A^ = (A^). 
Averaging was performed over not too faint stars, all epochs, and 
all parameters katR = /?opt, assuming that the mathematical ex- 
pectation of A^ does not depend on epoch, star light flux, and 
axis. This yielded the representative estimate 

A^ = 2(0.5V?' + 0.5^^' - 02)/ Yj^e'jDe). (41) 

i,e i,e 

For a complete set of FORS2 observations (any seeing con- 
ditions), and images in the central r < I'area, we obtained 
A - 160 //as. This value exceeds by much the standard deviation 
for epoch average residuals of bright stars, which is typically 
about 50-70 //as, and therefore such systematic eiTors strongly 
affect the histogram's shape (FigfT2b). It was found that large A 
values are associated with images of abnormally small FWHM, 
in particular FORS2 images in Dec 2006 with exceptionally 
good seeing of 0.3-0.5". Elimination of images with seeing be- 
low 0.47"(3.7 px) and the subsequent rejection of bad images 
with FWHM> 0.78"(6.2 px) significantly improved A. Such fil- 
tration produced better solution with essentially more compact 
histograms of epoch-normalized residuals (FigfTSb) fitted with 
a Gaussian of only 0.72 cr width parameter. This is insignifi- 
cantly smaller than the expected cr-width which, due to a change 
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Fig. 12. Histograms of Vg (epoch residuals Ve normalized to 
unity to exclude dependence on star brightness) in r = 1 'central 
zone (steps), Gaussian approximation (solid curves), and theo- 
retical distribution in the case of zero systematic errors (dashed 
Gaussians), for: a) - FORSl; b) - FORS2 complete set of im- 
ages, with seeing varying from 0.38"to 0.84"; c) - FORS2 subset 
of images with 0.47-0.78"seeing; d) - merged FORS2 residuals 
{V}e- Upper X-axis scale refers to the expected distribution of 
non-normalized epoch residuals for brightest stars of B=18 
mag for FORSl or R=16.5 mag for F0RS2. 



in N~^(i) caused by the above filtering, decreased from 0.78 to 
0.75. The described filtering of 20% of the FORS2 images, was 
applied to a few FORS 1 images whose seeing was almost always 
within the limits adopted. At the end of this Section we discuss 
the probable relation between image size and systematic errors. 

The frequency distribution of epoch average residuals {V)e 
based on merged residuals (V) (FigfTStl) does not differ from 
that built for V, (FiglH};). 

The estimates of A for stars in CCD central circular areas of 
r - 40", 1', and 1.5' are given in Table|4l 



Table 4. Systematic component A [//as] in residuals and in 
{V)e within central CCD zones of r radii 





FORSl 


FORS2 


r 


inV, in<V>, 


inV, in<V), 


40" 


53±21 25±19 


65±22 64±19 


1' 


61±9 69±10 


68±9 68±10 


1.5' 


95±10 79±12 


85+6 89±6 



According to Table [H the characteristics of the systematic 
error A for both cameras are similar and show a slight increase in 
the direction from the center to the periphery of a frame. Thus, 
while at the periphery A is near to 100 juas, at r < 40"it does 
not exceed 50-60 //as. Due to statistical limitations we not can 
estimate A at the center (where the target is usually placed), but 
considering the tendency observed we predict it could be about 
25 /(as at r < 10", as expected from the following discussion. 
The systematic components in {V}e and Ve are approximately 
equal. 

An important piece of information on the global distribu- 
tion of systematic errors over the CCD plane was derived apply- 
ing a low -pass Gaussian spatial filter to the epoch residuals Ve- 
The resulting low-frequency component A(x, y) in for each 
camera, some epochs, and reduction versions is shown in con- 
tour plots Fig ll3ll4l with isolines drawn with 25 /(as increments. 
All graphs refer to the residuals on the jc-axis. A{x,y) function 
change is rather complicated and has several extremums. This 
behavior, of course, not can be approximated by polynomials 
with basic functions / since this dependence is excluded in the 
course of the reduction procedure. It is characteristic that large 
systematic errors reaching in most cases 100-200 /(as tend to 
concentrate at the periphery. At the field center, where the target 
is usually placed, A{x,y) functions varies rather smoothly and 
often fall to < +25 //as. No graph was found to have extremum 
at the center of the frame. 
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Fig. 13. Global pattern of FORS 1 systematic error distribution 
over the CCD surface: a) - in normal points Ve and b) - in merged 
residuals {V} for a Dec 2002 epoch. Isolines are plotted every 
25 //as; zero level is shown by a solid line. Residuals Ve were 
computed with k - 10 at optimal size of reference frames 7?opt = 
1.5'. 



In the case of FORS 1 , systematic error plots are given for 
middle 2002 epoch at which the maximum fluctuations of A(x, y) 
are detected. FigfTSh. plotted for Ve computed with k = 10 at 
R = ^opt and FigfTJb for merged (V) residuals have few simi- 
lar structures above ±50 //as. At the center, errors are negligibly 
small. Although we discuss here only a particular case of the 
reduction with k - 10 or of merged (V), our comments (both 
for FORSl and FORS2) are valid also for computations with 
other k due to the high correlation of results obtained at /? = /?opt 
(Sect.4.3). 

FigfT4b shows the systematic pattern after elimination of im- 
ages with abnormally good and bad seeing and illustrates the 
decrease of systematic errors compared to the use of a complete 
set of images (FigfT4k). Few extrema have vanished and most 
contrast details are smoothed. These graphs correspond to the 
middle of the January 2007 epoch with the largest variance of 
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Fig. 14. Distribution over the CCD surface of a systematic com- 
ponent in FORS2 epoch residuals in Jan 2007: a) - in com- 
puted with A: = 10 and /?opt = 1.5'with all images; b) - the same, 
for a subset of images in a 0.47-0. 78"seeing range; c) - the same, 
for a full set of images and modified reduction model that takes 
into account dependence of positions on seeing (Sect.6.2); d) - 
for merged residuals (V); e) - in Ye at good seeing (Nov 2006); f) 
- in Ye (Jan 2007) for k-d and 7?opt - 0.5'. Isolines are plotted 
every 25 yuas; zero level is shown by a solid line. See explana- 
tions in the text. 



monthly residuals (FigfTTT). which was chosen to show the worst 
case of A{x,y). In November 2006 (first epoch) these eiTors are 
much better (FigfT4b). The error structure for the merged (V) 
residuals (FigfT4tl) in general is like that for Ye ( FigfT4b) at 
the same epoch. Figfl4lF refers to Ye computed at A = 6 and 
i?opt = 0.5'. The distribution shown here is similar to that for 
the reduction with A: = 10 (FigfT4b) and for (V) with ±50 yuas 
systematic errors at the center, but much larger peaks at frame 
boundaries. 

The dependence of A(x, y) on x, y could be a reason for ex- 
cluding systematic errors in a secondary iteration; however we 
had an insufficient number of reference stars. Alternatively, Ae 
can be treated as additional components of each star model pa- 
rameter ^ to be computed with other parameters. This however 
is not useful for planet search or microlensing applications due 
to the complete zeroing of useful astro me trie signals. 

Considering that a change of Ae between two adjacent 
epochs at the frame center is approximately 50 yuas (100-200 yuas 
within the entire FoV) for both cameras, we can estimate the 
stability of the FORS astrometric system over short time scales. 
Assuming that a given change occurs at about a month spacing. 



the daily rate of systematic change is 1-2 yuas (3-6 yuas for the 
whole FoV). Our previous study of FORS 1 errors (Lazorenko et 
al. 2007) have shown that a difference in the systematic compo- 
nent over a time scale of four days is either undetectably small 
or 30 yuas atmost. This is about 8 //as a day change in system- 
atic eiTors in the whole FoV, which is in accord with the cuiTent 
estimate. 

Recall that the amplitude of systematic errors of 50 yuas we 
referred to corresponds to poor observing conditions; the esti- 
mate of 25 //as is more relevant for normal conditions and targets 
at frame center. 

6.2. Instrumental background of systematic errors 

Based on the discussion in the previous subsection, we conclude 
that the characteristics of the systematic component A^. in magni- 
tude and in spatial behavior are identical for FORS 1 and FORS2. 
This is the second identity of these cameras derived based on ob- 
servations at rather different time scales, CCD type, and photo- 
metric bands (the first is a dependence of image motion variance 
on R), and is evidently due to identical optical design of the cam- 
eras. 

Systematic errors are probably generated during image cen- 
troiding due to the highly complicated star profiles, actually in- 
definite at the high accuracies at which we work. Consider that a 
typical eiTor of 100 yuas is only 2 x 10"^ fraction of the FWHM. 
Although FORS images are appropriate, the definition of the 
'image photocenter' for star profiles distorted by variable optical 
aberrations becomes uncertain. Therefore we use a centroiding 
procedure (Lazorenko 2006) specialized for finding the weighted 
photocenter of the image which is more stable to image defor- 
mations than the "profile center". 

At zero or constant geometric distortion, the position of the 
measured image centroid is subject only to random eiTors caused 
primarily by Poisson noise in the number of photons and by at- 
mospheric image motion. Systematic and extra random compo- 
nents in position appear in the following cases of instability: 

- Small deformation of a star profile (change of asymme- 
try) in time due to a slow change of optical abeiTations. This 
results in image photocenter shift proportional to the gradient 
of aberrations at the point. These shifts, correlated in x, y space, 
are detected (after a certain filtration introduced by astrometric 
reduction) as systematic residuals A(x,y); their change in time 
produces A component. The effect is lowest at frame center were 
images are most symmetric due to optimal optical performance. 

- A random change of image-to-image atmospheric PSF 
and seeing (shape and size). This affects image profiles both in 
size and shape and therefore shifts photocenters by an amount 
proportional to seeing fluctuation and to the measure of image 
asymmetry at the point, representing modulation of the opti- 
cal aberration field by a random signal. The observed effect is 
random in time and correlated in space, and thus mimics atmo- 
spheric image motion. In Sect.5 we classified it as "instrumen- 
tal" image motion with an amplitude exceeding that of atmo- 
spheric image motion and not dependent on exposure time. 

The combined effect of image asymmetry and seeing vari- 
ation is observed in an exaggerated form for stars with nearby 
companions, the light of which causes image asymmetry. It 
is easily detected as a linear dependence of Yi,n on FWHM, 
which in most bad seeing conditions produces enormous system- 
atic deviations exceeding by much the random errors and being 
grounds for discarding these measurements (Sect. 2). 

Systematic deviations of the average seeing at a single epoch 
from its average for a full set of data may produce systematic 
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bias of positions. To see whether this exists in our data, we per- 
formed reduction with an expanded set of parameters ^ adding 
two extra terms describing the linear dependence of a star's x,y 
position on seeing. With a full set (all seeing) of images, a new 
value of the systematic error A in an r < I'area decreased from 
its former 160 //as to 90 juas, and a much smoother shape of 
A{x,y) (FigfT3t) was obtained in comparison to the initial pat- 
tern (FigfTSh). The improvement, however, is seen primarily at 
the periphery where geometric distortions are largest. Thus vari- 
ations of seeing actually lead to systematic errors in epoch po- 
sitions. We however consider the use of the expanded reduction 
model insufficiently validated, since it reduces the useful astro- 
metric signal and, for centrally placed targets, does not leads to 
a significant improvement in the accuracy. 

Thus, variations in seeing are a source of systematic errors; 
for precise astrometry these variations should be limited. 

7. Allan precision 

Another characterization of the VLT long-term astrometric sta- 
bility is based on the computation of the Allan deviation of resid- 
uals y,„. This quantity is normally used as a powerful indicator 
of systematic errors in observations and corresponds to the as- 
trometric precision of a time series of n images. FigfTSl presents 
plots of this variable as a function of time lag (expressed as the 
number of images «) between subsamples of residuals V,,,. The 
data in Fig[T5]are the average of each star Allan deviation taken 
over all, except the most peripheral, stars. Before averaging, we 
normalized the individual Allan deviations to their values at zero 
time lag so as to compensate for the highly varying amplitude of 
this variable for stars of different brightness. The right vertical 
axis of Fig {15] is the scale for the normalized Allan deviation 
computed as described and is valid for stars of any magnitude. 
The left axis is the scale used to find the Allan deviation for 
bright stars (B=18 mag at FORSl and /;=16.5 mag at F0RS2) 
with e = 230 juas, 6 - 150 juas and assuming its location at the 
frame center, which ensures small Ajf = 150 //as. 

For computations, we used V,,, residuals obtained with each 
reduction parameter k - 6... 16 and R - /?opt- Results for dif- 
ferent k are very similar due to the high degree of correlation 
between these data sets (dashed curves in Fig[T5]l; for FORS2 
they actually run into a single line corresponding to the Allan 
deviation of residuals (V). Besides, due to the normalizing pro- 
cedure, computations based on a subset of either all or bright 
stars only produced similar estimates. The Allan deviation was 
found to follow „-0-522±o.oo2 (FORSl) and „-0-506±o.oo4 (FORS2) 

power laws which are near to that expected for the average of 
a random variable. A simple comparison with an n"'^^ law is 
of course incorrect since least square fit residuals are correlated 
and have a non-diagonal covariance matrix B. To compare our 
results with those expected with zero systematic errors, we per- 
formed a numerical simulation of observations introducing un- 
corrected random noise in model measurements. The obtained 
dependence (solid curves in FigfTSll follow a power law with a 
slope n""^^^=^'"*'^ for both cameras, which is near to that ob- 
tained from observations. 

The difference in the observed and expected plots is clearly 
seen for the FORS2 camera starting from n > 15. This diver- 
gence is related to systematic errors discussed in Sect.6. Due to 
the small magnitude of the errors, they cause only a 5-10% in- 
crease in the astrometric error in comparison to that expected in 
the absence of systematic errors. 

Above, in the case of FORS2, we used images with a stan- 
dard FWHM to avoid any degradation of precision caused by im- 




Time lag (number of frames n) 

Fig. 15. Normalized Allan deviation (right axis) in positional 
residuals and astrometric precision for bright targets (left axis) 
expected from a series of n images with standard seeing only. 
Estimates are based on: - residuals V obtained with k - 6... 16 
(dashed lines which for FORS2 actually run into a single line); - 
merged residuals (V) (open circles); - numerical simulation as- 
suming zero systematic errors (thick line); - FORS2 residuals 
with no restriction on seeing (black dots). 

ages with abnormal seeing. The Allan deviations with all avail- 
able measurements used (black dots in FigfTSll are seen to have 
an excess of about 30% for n > 15. 




R magnitude 

Fig. 16. Astrometric precision (Allan deviation) for a series of 
30 images as a function of magnitude derived from five-month 
FORS2 observations (open circles); the same expected in the 
case of a very dense reference frame (Arf — 0) and moderate 
6 - 150 /las image motion (solid curve). 

Given a series of 30 images (0.5-1 hours of telescope time), 
the precision of FORS 1/2 astrometry (Fig[T5]l is about 50 //as. At 
this fixed number of images, Fig[T6](open circles) shows the de- 
pendence of the precision on star brightness computed for each 
star as the individual Allan deviation. Data refer to FORS2 for 
which sufficiently long observation series are available and for 
images obtained at normal seeing. Because precision very much 
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depends on the reference frame noise (especially for bright tar- 
gets), the plot shows only those stars for which Arf < 300 //as. 
For that reason, many of the brightest 16.5-17.0 mag stars were 
omitted due to their peripheral location in the frame. The ob- 
served dependence of precision on magnitude shows it to be bet- 
ter than 100 fias for R <19 mag targets and about 50 yuas for 
R 11 mag. 

The above estimate includes the error caused by reference 
frame and optical aberrations which depend on a particular den- 
sity of star distribution in the sky, stability of the optical system, 
and variations of seeing conditions at the period of observations. 
Precision improves with improved conditions. A case of interest 
is the precision expected at high reference star density (Ajf = 0) 
and image motion 5 = 150 juas, typical for a reference frame size 
of 0.7-1.5'. These estimates, shown by a solid line in FigfTSl 
prove the feasibility of 50 fias astrometric precision for brightest 
targets. 

8. Conclusion 

Astrometric quantities (residuals of positions, model parame- 
ters) derived from the processing of images are intrinsically rel- 
ative. They are computed in a certain system and relative to 
the frame of reference objects specific to the particular target. 
These peculiarities of differential reduction should be taken into 
account for the interpretation of the output data and in more 
complicated cases of handling inhomogeneous series of images, 
for instance obtained in different spectral bands or even with 
different cameras, which is a case expected for long-term pro- 
grammes. We demonstrated that a careful processing of B and R 
images in a common system does not degrade precision. 

The precision of astrometric imaging at VLT depends on sev- 
eral noise sources. Uncertainty of the image photocenter deter- 
mination e, of course, is the dominant component of the total 
eiTor Our data show that Eq.([T]) provides a correct estimate of 
this error, at least with an accuracy of +5%, for a wide range 
of light fluxes and seeing conditions, s depends not on the star 
magnitude but on the light flux collected in the star image. For 
that reason, the lowest s (equally, best astrometric precision) is 
expected for images at saturation level, which depends on the ex- 
posure, filter, pixel scale, and seeing. Therefore, having, for ex- 
ample, dependence s{R) of s on magnitude R, it is easy to apply 
it to other observations. For instance, in B band s{B) x E(R+Am) 
where Am is the difference of magnitudes of images with equal 
light flux in B and R bands. This scaling is illustrated in Fig|9] 
where s(B) and s(R) are seen to be the same dependences shifted 
by Am ^ 1.5 mag. 

Mitigation of atmospheric image motion at T > 70 s expo- 
sure does not present a problem due to its small amplitude in 
comparison to other random noise components. Quite unexpect- 
edly, however, we found that the measured image motion vari- 
ance is the same at the very different T - 70 and T x 600 s ex- 
posure. We consider this as a new type of random error caused 
by the combined effect of the telescope-related asymmetry of 
star profiles and of random changes in the atmospheric PSF. We 
have found that images with abnormally bad and good seeing 
are affected by large random and systematic errors and thus their 
use should be avoided in precision astrometry. 

Estimates of astrometric eiTors obtained in this study re- 
fer rather to precision caused by stochastic error components 
but not to accuracy. The final astrometric accuracy can de- 
grade compared to precision due to systematic errors incorpo- 
rated into the reference frame and then propagated by the reduc- 
tion model. Presently, our astrometric results cannot be cross- 



compared against independent datasets and analyses of common 
objects, which could provide direct estimates of the long-term 
accuracy. In Sect.6,7, however, we demonstrated that the mea- 
sured star displacements in space and time are fitted to nearly 
white noise residuals with a variance predicted by the model, 
therefore, an extra systematic component (if present) should 
closely trace both the parallactic and proper motion of each star. 
A too low probability of this scenario (considering large accu- 
mulated Dof) implies that the accuracy of the VLT relative par- 
allax and proper motion determination is comparable to the pre- 
cision. 

We have demonstrated that, with reference to systematic er- 
rors, astrometry at FORS 1/2 is accurate to 25 juas at five month 
and 6 year time intervals. Due to that fact, relative proper mo- 
tions and trigonometric parallaxes of stars in the test field were 
derived with a precision of 20 //as yr"' and 40 //as respectively 
for 17-19 mag stars. Thus distances to stars at 1 kpc can be mea- 
sured with a precision of 4%, providing a coiTection from rela- 
tive to absolute parallax is added. This level of accuracy and, 
especially, good long-term stability, has numerous astrometric 
applications including measurement of astrometric microlens- 
ing, planet detection and characterization by measuring reflex 
motion of the parent star, and kinematics of Galactic stellar pop- 
ulations. In the context of exoplanet searches, the use of the VLT 
to search for planets near brown dwarfs is very efficient. With a 2 
year observation programm, Saturn mass planets with orbit peri- 
ods longer than 1 year and Netpune mass planets with 2 year pe- 
riods are detectable for brown dwarfs located at 10-20 pc. These 
observations will make a real breakthrough in our understanding 
of planet and brown dwarf formation by probing a separation 
and mass range that is poorly suited to other techniques. 

CuiTently our method has several crucial limitations. It is 
applicable to moderately populated sky areas with low rates of 
blending but that are rich enough to provide a sufficient num- 
ber of reference stars. Also, targets brighter than 15-16 mag not 
can be measured due to saturation if the exposure time is not too 
short. Therefore it is not applicable to highly crowded sky ar- 
eas like those used for microlensing works and to bright, nearby 
solar-type stars which are objects of interest for planet searches. 
Most of these limitations originate from the star profile fitting 
and can be removed with use of a better technique for determina- 
tion of star photocenters. Precision astrometry is also problem- 
atic for telescopes with segmented primary mirrors that produce 
intrinsic complex and time-variable PSF. 

We emphasize that the results of this study do not refer to a 
specific telescope, they validate precision astrometry for a whole 
class of large ground-based imaging telescopes. Astrometry ben- 
efits highly from the use of large apertures, mitigating in this way 
both the principal image centroiding component of the total error 
and atmospheric image motion. For future 30 m telescopes, the 
astrometric precision is expected to be 10//as or better for image 
series of 0.5 hour duration. At this precision, we are challenged 
with a wide spectrum of problems, e.g. the ability to measure 
bright objects due to saturation and the predominance of sys- 
tematic errors caused by optical aberrations. 
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